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Abstract 

We examine the adjacency matrices of three-regular graphs rep- 
resenting one-face maps. Numerical studies reveal that the limiting 
eigenvalue statistics of these matrices are the same as those of much 
larger, and more widely studied classes from Random Matrix The- 
ory. We present an algorithm for generating matrices corresponding to 
maps of genus zero, and find the eigenvalue statistics in the genus zero 
case differ strikingly from those of higher genus. These results lead us 
to conjecture that the eigenvalue statistics depend on the rigidity of 
the underlying map, and the distribution of scaled eigenvalue spacings 
shifts from that of the Gaussian Orthogonal Ensemble to the exponen- 
tial distribution as the map size increases relative to the genus. 



1 Background 

In this paper we examine a special class of three-regular graphs. A graph 
G = {V, E} consists of a set of vertices V and a set of edges E, with each 
edge connecting two vertices. A graph is fc-regular if the number of edges 
incident to each vertex is exactly k. The structure of a graph is encoded in 
the associated adjacency matrix A = [uij], where Ojj is the number of edges 
between vertex i and vertex j. When we refer to the spectrum of a graph, 
we are referring to the eigenvalues of its adjacency matrix. 

Our special class of three-regular graphs have 2N vertices and can be de- 
composed into a cycle graph and a totally disjoint, one-regular "toothpick" 
graph. The cycle graph on 2N vertices is the two-regular graph in which 
vertex i is adjacent to vertices i + 1 and i — 1 mod 2N. In the standard 
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toothpick graph, vertex 1 is adjacent to vertex 2, vertex 3 is adjacent to 
vertex 4, and so on. The adjacency matrices of our graphs can be expressed 
as C + P^TP, where C is the adjacency matrix of the cycle graph, T is the 
adjacency matrix of the standard toothpick graph, and P is a permutation 
matrix. Three-regular graphs of this form can be used to represent one-face 
maps. 

Definition 1.1. A map (or dessin d'enfant) is an embedding of a connected, 
labeled graph X , into a compact oriented surface S, with the following prop- 
erties: 

1. The edges of the graph do not intersect (except at vertex points). 

2. The complement of X in S is a disjoint union of open cells which are 
all homeomorphic to the disk. These cells are called the faces of the 
map. 

The genus of a map is the genus of the underlying surface, not the graph- 
theoretic genus of the embedded graph X. 

Our three-regular graphs describe a glueing of the 2A^-gon's edges. If the 
toothpick portion of the three-regular graph contains an edge connecting 
vertex i and vertex j, then edge i and edge j are glued together in the cor- 
responding 2A^-gon. This glueing maps the 2A^-gon to a graph embedded 
in a Riemann surface of genus 5 > 0. Since the complement of this graph 
is homeomorphic to the unit disc, we say the map has "one face". Figure 
[T] shows a one-face map, its representation by a 2A^-gon with edge identifi- 
cations, and the associated three-regular graph. All one-face maps can be 
represented by a three-regular graph of the form C + P^TP. 
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Figure 1: a. The genus one, one-face map b. The corresponding 2A^-gon 
with edge identifications c. The associated three-regular graph of the form 
C + P^TP 
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2 Numerical studies of one-face maps 
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Figure 2: Empirical eigenvalue density for a random sample of 300 matrices 
having the form C + P^TP, plotted with the McKay density function. The 
corresponding graphs have 600 vertices each. 



We generated large samples of random matrices having the form C + P TP, 
using Matlab's random permutation generator. The limiting eigenvalue 
statistics of our random matrices are the same as those of much larger, 
and more widely studied classes. In particular, the limiting eigenvalue den- 
sity is described by the McKay density [T], and the distribution of scaled 
eigenvalue spacingqll from the bulk of the spectrum appears to be that of the 
Gaussian Orthogonal Ensemble (GOE). Figure [2] shows the empirical eigen- 
value density for a random sample of C + P^TP matrices. Figure [3] shows 
the scaled eigenvalue spacing distribution over the bulk of the spectrum for 
the same sample. 

The McKay density formula, 

= ( f« l-l < I (1) 

[ otherwise J 

describes the limiting eigenvalue density for sequences of of /c-regular graphs 
which look like trees at most vertices. 

Theorem 2.1 (McKay). Given a sequence Xi, X2, . . ■ of k -regular graphs 
for some k > 2, let \Xn\ denote the number of vertices in graph X^, and 

^Spacings are found by ordering the eigenvalues in ascending order and finding the 
distance between consecutive eigenvalues. These distances are then scaled to have mean 
one and integrated density one. 
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Figure 3: Distribution of scaled eigenvalue spacings over the bulk of the 
spectrum for a random sample of 300 matrices having the form C + P^TP, 
plotted with the GOE bulk scaling limit. The corresponding graphs have 
600 vertices each. 



Wr{Xn) the number of closed walks 
and 

Wr{Xn) 
n^oo |A„| 

then the limiting probability density 
Xn is fkix). 



of length r. If — > co as n ^ oo, 

= Vr > 3 (2) 
function describing the eigenvalues of 



McKay proved that for a deterministic sequence of A:-regular graphs, if the 
limiting eigenvalue density is not described by McKay's density formula, 
then the ratio of the number of closed walks of length r to the number of 
vertices must not approach zero, for some value of r. 

The Gaussian Orthogonal Ensemble is the set oi N x N matrices whose 
components are independent (up to the symmetry requirement) random 
variables of mean zero and standard deviation one. Numerical studies have 
revealed that the distributions describing the eigenvalue statistics of the 
Gaussian Orthogonal Ensemble (GOE), describe the eigenvalue statistics of 
most real symmetric random matrices having independent entries. How- 
ever, the entries of our random matrices are highly dependent. Since these 
matrices represent three-regular graphs, the entries of each row and each 
column must sum to three. Furthermore, the cycle subgraph ensures that 
two of the non-zero entries of each row and each column come from the first 
sub / super-diagonals of the matrix. 

Given our random matrices describe one-face maps, a natural question is 
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whether the eigenvalue statistics depend on the genus of the underlying 
map. To investigate this question we must generate matrices corresponding 
to random one- face maps having N edges and genus g. 



3 Constructing genus zero maps 

The number of genus g one-face maps having N edges in the embedded 
graph is given by [21 Theorem 2] 



eJN) = — X the coefficient of x^^ in ( — ^/'^ , ^ 

^ {N + 1)\{N - 2gy. Vtanhx/2y 

(3) 

The vast majority of one- face maps having edges are of genus slightly 
less than [N/2\ . For a given N, random selections of one-face maps will 
only sample genera near this upper bound. For large A^, the distribution of 
genera makes it infeasible to generate a random sample of low genus (relative 
to [N/2\) by generating random one-face maps and sorting out the maps 
with desired genus. We need an algorithm which will generate random one- 
face maps of a specified genus. In the genus zero case, the bijection between 
one-face maps and pair partitions provides such an algorithm. 

As previously described, the important feature in the representation of one- 
face maps by three-regular graphs of the form C + P^TP is the toothpick 
subgraph, P^TP. This portion of the graph defines the edge glueings of 
the 2A'^ polygonal representation. This information can be encoded in a pair 
partition diagram. If vertex i and vertex j are adjacent in the toothpick 
subgraph, then node i and node j are paired in the associated pair partition. 
There is a bijection between the set of non-crossing pair partitions and genus 
zero one-face maps. Using this bijection, we have developed an algorithm to 
generate the adjacency matrices of three-regular graphs representing random 
genus zero one-face maps. The generation of one-face maps having specified 
genus g > and A^ edges remains an open problem. 

Our algorithm generates a random non-crossing pair partition and then uses 
the node identifications of this pair partition to assign adjacencies in the 
toothpick subgraph of the associated three-regular graph. Non-crossing pair 
partitions are formed recursively. The node paired with node one is ran- 
domly chosen according to the probability mass function described below. 
This pairing partitions the remaining nodes into two sets, one within the 
pairing and one external to it. The first node of each of these sets is then 
paired using the same algorithm. This process is repeated until all nodes 
have been paired. 
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Figure 4: Three-regular graph representing a genus zero one-face map and 
the associated non-crossing pair partition 



The number of non-crossing pair partition on 2N nodes is given by the 
Catalan number, 

„ . 1 f2N- 
Un ■- 



N+1\N 
which can be defined recursively as 

_ (4iV - 2) • Cn-1 

- — wn — • 

Let 2m be the random variable representing the node paired with node one 
in a non-crossing pair partition. Since the pair partition is non-crossing, 
odd nodes must be paired with even nodes. The number of non-crossing 
pair partitions on 2N nodes in which the first node is paired with node 2m 
is equal to the number of non-crossing pair partitions on 2m — 2 elements 
times the number of non-crossing pair partitions on 2N — 2m elements, 
Cm-i ■ C^-m- Thus, the probability mass function {pmf) of m is a function 
of both m and N, and is given by 

pmf{m,N) = ^"i-i CN-m ^ 

To avoid computer arithmetic errors, the pmf{m,N) for m greater than 
five and less than is calculated using the recursive formula derived from 
equation 

Our algorithm was implemented using Matlab. With each recursive call 
to our program, we assign a pairing to the first element of some subset of 
the entire pair partition. These subsets are determined by the previously 
assigned pairs. For the example shown in Figure [5l after the first pairing is 



6 



assigned, the partition is divided into two subsets of lengths 2mo — 2 and 
2N — 2mo- The first element of each of these subsets is paired during the 
next round of calls. Since the pmf used to assign these pairings is based on 
the subset lengths, it is different for each recursive call. 




Figure 5: a. The partially completed pair partition after one call to our 
algorithm. The mo is randomly chosen according to the pmf ^ and node 1 
is paired with node 2mo. b. The partially completed pair partition after a 
second round of calls. Elements mi and m2 are randomly selected according 
to the probability mass functions pmf{m,mQ — 1) and pmf{m,N — mo), 
respectively. 



3.1 Proof that genus zero one- face maps are chosen with 
equal probabihty 

As previously discussed, randomly choosing a genus zero one-face map hav- 
ing edges is equivalent to randomly choosing a non-crossing pair partition 
on 2N nodes. The number of non-crossing pair partitions on 2N nodes is 
Cn. Thus, we want our algorithm to construct non-crossing pair partitions 
on 2A^ nodes in such a way that every partition has a ^ chance of being 
created. 

Induction Case: There is only one way to partition 2 nodes. Node 1 is paired 
with node 2 with probability 1, and the resulting non-crossing pair partition 
is chosen with probability 1 = = j. 

Suppose that for all k less than N, the algorithm randomly constructs non- 
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crossing pair partitions on 2k nodes in such a way that each possible partition 
has the equal probability of being created, Let X be an arbitrary non- 
crossing pair partition of 2N elements, and let node 2m be the element paired 
with node one in X. The probability our algorithm constructs partition X 
is equal to the probability of node 1 being paired to node 2m, times the 
probability of the particular non-crossing pair partition / on the 2m — 2 
nodes internal to the pairing 1:2m, times the probability of the particular 
non-crossing pair partition E on the 2N — 2m nodes external to the pairing 
1:2m. By our supposition, our algorithm creates the internal pair partition 
/ with probability ^ ^ , and the external pair partition E with probability 
— . Thus, the probability of an arbitrary non-crossing pair partition X 
being created is equal to 

P{X) = pmf{m,N) ^ ^ 



Cm-l C]Sf-m 

Cn Cm-i Cn- 

1 

Cn 



(7) 



By induction, our algorithm randomly constructs non-crossing pair parti- 
tions on nodes with equal probability 



4 Results 

We can examine the eigenvalue statistics of one-face maps having specified 
genus near [N/2\ by forming large samples of random one-face maps having 
A'' edges and sorting by genus. Figures [6] and [7| show the eigenvalue density 
and spacing distribution for a random sample of genus 147 one-face maps 
generated in this way. Examining the eigenvalue statistics of one-face maps 
having genus near [N/2\ , we find that for large the eigenvalue density ap- 
pears to be McKay, and the scaled spacing density appears to be that of the 
Gaussian Orthogonal Ensemble. These results are unsurprising given that 
samples of random one-face maps for large N exhibit the same properties 
and are comprised primarily of maps having genera near \_N/ 2j . 

Using the previously described algorithm, we generated samples of random 
genus zero one-face maps for large N, and examined the eigenvalue statis- 
tics. The results are strikingly different from those found for one- face maps 
having genus near [N/2\ . The distribution of eigenvalues for a random 
sample of genus zero one-face maps is plotted in Figure El In Figure [S] the 
eigenvalue density curves for various samples of genus zero one-face maps 
are plotted on the same axis. Although the number of embedded edges A'^ is 
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Figure 6: a. The eigenvalue density of three-regular graphs representing 
genus 147 one-face maps. The sample was created by generating random 
three regular graphs of the form C + P^TP corresponding to one-face maps 
having N = 300 edges. The graphs were then sorted according to the genus 
of the associated map. Thirty-six genus 147 maps were included in this 
sample, b. The eigenvalue density of three-regular graphs representing genus 
zero one-face maps. The sample contains 100 graphs, each representing a 
genus zero one-face map having = 1000 edges. 



different for each sample, the eigenvalue densities are nearly identical. While 
this distribution is the same for all large values of A^, it does not appear to 
be described by any established density formula. Examining the eigenval- 
ues themselves, we find that the peaks observed in the density curve near 
0, ±0.5, ±1.8 and ±2.3 are not the result of repeated eigenvalues, but rather 
a decrease in spacing between consecutive eigenvalues. The non-crossing 
structure of the toothpick portion of the three-regular graphical represen- 
tation of genus zero one-face maps ensures that the graphs are bipartite. 
The symmetry in the distribution of eigenvalues about A = is due to this 
bipartite structure. 

The difference in the distribution of eigenvalues for the genus zero and genus 
~ [N/2\ cases is also seen in the corresponding mean j'*'' spacing curves. 
Figure [9] shows the mean j^^ eigenvalue spacing as a function of j for a 
sample of three-regular graphs representing genus zero one-face maps, and 
a sample of three-regular graphs representing one-face maps of genera near 

LfJ. 

It is clear from these results that, in the genus zero case, the distribution of 
eigenvalues is not described by the McKay density formula. In fact, three- 
regular graphs representing genus zero one face maps violate the criteria 
of McKay's theorem (12. ip . These graphs are always planar, i.e. they can 
be drawn in the plane in such a way that the edges of the graph do not 
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Figure 7: a. The scaled eigenvalue spacing distribution of three-regular 
graphs representing genus 147 one-face maps plotted with the GOE bulk 
scaling limit. The sample was created by generating random three regu- 
lar graphs of the form C + P^TP corresponding to one-face maps having 
N = 300 edges. The graphs were then sorted according to the genus of the 
associated map. One-hundred and nineteen genus 147 maps were included in 
this sample, b. The scaled eigenvalue spacing distribution of three-regular 
graphs representing genus zero one-face maps, plotted with the exponential 
probability density function. The sample contains 100 graphs, each repre- 
senting a genus zero one- face map having N = 800 edges. 
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Figure 8: The eigenvalue density curves for random samples of genus zero 
one-face maps having N = 50, = 200, and = 400 edges. Each sample 
contains 500 genus zero one-face maps. 
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Figure 9: a. The mean j*^ eigenvalue spacing as a function of j for a 
sample of 100 three-regular graphs representing one-face maps of genera 
near g = 250. b. The mean j'*'' eigenvalue spacing as a function of j for a 
sample of 100 three-regular graphs representing genus zero one-face maps. 
Both samples contain graphs having 1000 vertices. 



intersect except at vertex points. A planar drawing of a graph divides the 
set of points of the plane not lying on the graph into regions called faces. 
Figure [10] shows the three-regular graph representing a genus zero one- face 
map with its faces labeled. One of the faces will be unbounded. This face 
is called the infinite face. All other faces are bounded by the edges of the 
graph, and are called internal faces. 

For genus zero one-face maps, the embedded graph is a planar tree, and each 
internal face of the associated three-regular graph corresponds to a vertex 
of this embedded tree. The degree of this embedded vertex is equal to ^ the 
face boundary length. Thus, the distribution of vertex degrees in the genus 
zero embedded trees determines the distribution of face boundary lengths in 
the associated three-regular graph. Furthermore, since each face boundary 
forms a cycle in the graph, the number of vertices of degree d in the tree 
provides a lower bound for the number of 2d cycles in the three-regular 
graph. 

Let X]sf denote the set of three-regular graphs having 2N vertices which 
represent genus zero one-face maps. Drmota and Gittenberger [3| proved 
that the number of degree k vertices in a random planar tree having A'^ 
edges is asymptotically Thus, the expected number of 2fe-cycles in a 
random element of as A^ ^ oo is 

A^ 

lim E(c2fc) > 

Since the number of closed walks of length 2k is greater than or equal to the 
number of cycles of length 2fc, the set of three-regular graphs representing 
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Figure 10: The three-regular graphical representation of a genus zero on-face 
map with its faces labeled. 



genus zero one- face maps contains a subset {Xi,X2, ■ ■ ■} of nonzero measure 
for which 



This violates the criteria of McKay's Theorem. 

The distribution of scaled eigenvalue spacings for a sample of random genus 
zero one-face maps having N = 800 edges is plotted in Figure [71 The 
scaled spacings in the genus zero case are not described by the GOE density 
formula which was shown numerically to describe the bulk scaling limit for 
genera near [N/2\, but rather appear to be described by the exponential 
probability density function. 

5 Conclusion and conjecture 

The numerical results suggest that the eigenvalue statistics depend on the 
rigidity of the associated three-regular graph. For one-face maps having 
small genus relative to A^, a great deal of structure is imposed on the corre- 
sponding three-regular graph. We believe the striking difference between the 
eigenvalue statistics of genus zero maps and genus ~ [X'/2\ is a result of this 
imposed structure. For any fixed genus, as the size of the map increases (i.e. 
the number of vertices in the associated three-regular graph grows), more 
structure will be imposed on the associated three-regular graph. This leads 
to the conjecture that for fixed genus g, in the limit as ^ oo, the scaled 
spacing distribution will be exponential. As we have seen, for N ~ 2g, the 
scaled spacing distribution appears to be GOE. Thus the distribution must 
undergo a shift from GOE to exponential as N increases. Similar shifts in 
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probability density functions have been observed in biological networks and 
quantum chaotic systems [U [5] . 
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